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We calculate the charge distributions of d*(2380). Two different interpretations of the d* are 
considered for a comparison. One is a compact explanation with coupled AA + CC two-channel 
approximation in the chiral constituent quark model. Another is a resonance state of Dir. The 
remarkable differences of the charge distributions in the two pictures are shown and it is expected 
that the future experiments may provide a clear test for the different theoretical interpretations. 
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I. INTRODUCTION 


d* (2380) is a new resonance recently observed by CELSIUS/WASA and WASA@COSY Collaborations [I], [2j. It 
was found in the analysis of double pionic fusion channels pn — dr°r° and pn — drn when the ABC effect 
and the analyzing power A, of the neutron-proton scattering data were studied. It is argued that the observed 
structure cannot be simply understood by either the intermediate Roper excitation contribution or by the t-channel 
AA process, Refs. [1| [2] proposed an assumption of existing a d* resonance whose quantum number, mass, and 
width are I(J?) = 0(3+), M ~ 2370 MeV and I’ = 70 MeV (see also their recent paper [4], the averaged mass and 
width are M ~ 2375 MeV and I = 75 MeV, respectively). Since the baryon number of d* is 2, it would be treated 
as a dibaryon, and could be explained by either ”an exotic compact particle” or ”a hadronic molecule state” 
Moreover, since the observed mass of the d* is about 80 MeV below the AA threshold and about 70 MeV above the 
ATN threshold, the threshold (or cusp) effect may not be so significant as that in the XYZ particles (see the review 
of XYZ particles lel for example). Thus, understanding the internal structure of d* would be of great interest. 


The existence of such a non-trivial six-quark configuration with J(JP) = 0(3*) (called d* lately) has triggered 
a great attention and has intensively been studied in the literature even before the COSY’s discovery . In 
fact, after the experimental observation of the d*, there are mainly three types of model explanation for its nature 
in the market. The first one proposed a AA resonance structure and obtained a binding energy of about 71 
MeV (namely Mas = 2393 MeV) and a width of about 150 MeV. The second one announced a broad resonant 
structure of Dim for d*, whose resonant pole is around (2363 + 20) +i (65 + 17) MeV. The third one [15], following 
our previous prediction (ay, suggested a dominant hexaquark structure for d*, with a mass of about 2380 — 2414 MeV 
and a width of about 71 MeV, respectively Hig. From the results of those models, one has three observations: 
(1) All proposed models can reproduce a right mass of d*; (2) Only can last two models provide a total width for 
d* which compatible with the observed data; (3) only can the last model give all partial decay widths whose values 
agree with the observed data quite well. Even more, in terms of the third model (hexaquark dominant model), 
the predicted decay width of the single pion decay mode of d* — NNr is about ”1 MeV” which is much smaller 
than the double-pion decay widths fe} Here, we would particularly mention that this value is much smaller than 
that with the second scenario of Dior, but agrees with the experimental observation very recently [g]. All the 
outcomes from the third scenario support the idea that d* is probably a compact six-quark dominated exotic state 
due to its large CC component by Bashkanov, Brodsky and Clement [20}. Therefore, the model decay width of the 
d* — NWNrz process could be one of the way to justify the structure of d*, namely the width would be sensitive to 
the structure of d*. Although we have this weapon in hand, we are now still facing the problem that there is any ad- 
ditional way to justify different structure models. The general review on the dibaryon studies can be found in Ref. 2a}. 


It is known that the electromagnetic probe is one of the most useful tools to test the internal structure of a 
complicated system. For example, the electromagnetic form factors of the nucleon show the charge and magnetron 
distributions of a nucleon. The slops of the charge and magnetic distributions at original gives the charge and 
magnetic radii of the system. The precise measurement of the proton charge radius provides a criteria for different 
model calculations. For the spin-1 particle, like a deuteron or a p- meson, the charge, magnetic and quadrupole form 
factors tell the intrinsic structures as well, like charge and magnetron distributions and the quadrupole deformation 
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of the system. Therefore, the form factors of d*, for instance the charge distribution, might also be a physical 
quantity for discriminating the structure of d*. In this work, the charge distribution of the new resonance d* (2380) 
will be discussed with two different theoretical scenarios for comparison. One is a system with a single AA structure 
or a coupled AA + CC structure in our chiral constituent quark model. The other one is a resonant system of Dior. 
Moreover, the d*(2380) is a spin-3 particle, it has 25 + 1 = 7 form factors. A detailed discussion of all the seven 
form factors is beyond the scope of this work and will be given elsewhere. Here we only concentrate on its charge 
distribution and the charge radius of the d*(2380) in the two scenarios. 


This paper is organized as follows. In Sect. II, a brief discussion about the electromagnetic form factors of particles 
with spin-1/2, spin-1, and spin-3 will be shown. An explicit calculation of the charge distribution of d* with two 
scenarios is given in Sect. III. Sect. IV is devoted to a short summary. 


II. ELECTROMAGNETIC FORM FACTORS 


The study of the electromagnetic form factors of the nucleon (spin-1/2) is of great interest because it can tell us 
the information about the charge and magnetron distributions of a nucleon. In the one-photon approximation, the 
electromagnetic current of a nucleon is 


< NW) | Jk | NO) >= Ow) [Ph + i POVO), (1) 


where My is nucleon mass, q = p’ — p is momentum transfer, Q? = —q?, and Fı (Q?) and F(Q?) are the Dirac and 
Pauli form factors, respectively. These two form factors relate to the electric and magnetic form factors 


Gr(Q?) = F(Q’) — nFo(Q?), = Gur (Q?) = F(Q?) + P(Q’), (2) 


with n = Q?/4M3,. The normalization conditions of the two form factors for the proton and neutron are G% (0) = 1, 
GEO) =0, Gu = 2. 79, and G% = —1.91, respectively. 


In the Breit frame, we have q” = (0,@), p? = p? = M}, p=—p = —49, and po = ph) = E = My(1+7n)'/?. Then, 
we obtain 
< N(G/2) | Jn | N(-@/2) > = (1+ nm PxqixsGe(Q”) (3) 
E > E LOX 
< NGD | Jn ING) > = 0t) h Eau. 


Clearly, the electric form factor G'g (Q?) is directly related to the matrix element of J9. 


For a spin-1 particle, like deuteron, it contains three form factors. In the one-photon approximation, the electro- 
magnetic current is 


Jip (p, p) =en (p NSE gee (P ) (4) 


where e and ¢' stand for the polarization vectors of the incoming and outgoing deuterons, i and k are the polarizations 
of the two deuterons, and 
gu Qa Qg 
2 
mp 


ha = —[C1(Q)gas - 6310P) ZEE] P” — Ga(Q?)(Qagh — Qag) , (5) 


where P = p' +p. The three form factors G1,2,3(Q°) relate to the charge Gc(Q?), magnetic Gm (Q?) and quadrupole 
form factors Go (Q?) as 


Go(Q) = G(Q?) + ;mvGl@),  Gu(Q?) = GQ?) , (6) 
Gal@?) = GilQ) — G(Q?) + (1+ nv )GaQ@) , (7) 


with np = Q? /4m3 and Mp is the deuteron mass. The charge, magnetic and quadrupole form factors are normalized 
to Go(0) = 1, Gu (0) = 2m4 = 1.714, and Ge(Q?) = mpQa = 25.83, respectively. In the Breit frame, one may 
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also see that the charge form factor of the deuteron Go(Q?) can be obtained by direct calculating the matrix element 
$y <P,A| J? | prA>. 


For the d*(2380) particle, since its spin is 3, it has 2s+1=7 form factors. Its field can be expressed as €agy a rank-3 
tensor which is traceless. Clearly, €aag = 0, €agy = EBay, and p“€agy = 0. In the one-photon exchange approximation, 
the general form of the electromagnetic current of a 37 particle is 


7 a’ B’ y. a 
weiter (PMi gry" apy BY (p) (8) 


and the matrix element 


Merpry aby = |Gi(Q*)P# [Gara (98897 + G8'494'8) + permutations] (9) 


+G2(Q 
+G (Q 
+G4(Q 
pa 


+Ge(Q”) 


2)PH | garda [9889y + 99'94'8] + permutations] /(2M?) 


Pype [aw gaqoraogya + permutations] /(4M*) 
P” qaqaqa qedy'qy/(8M°) 
| (9h = gia") C T 9993") + permutations] 


a 


(otqa = gga ) (araso + Gy! Gy 98'B + WB! Gy Gy'B + ay asgye ) 


+ permutations] /(2M?) 


+G7(Q’) (ota = ghia") 48'98 Vy! Vy + permutations] /(4M | 


9479,4,9),0, 


invariant condition 
quM ar) aby =0 (10) 


is fulfilled as well as the time-reversal invariance. The combinations of the above seven form factors Gi2,...7(Q?) can 
give the physical form factors of d* such as the charge, magnetic, quadrupole as well as other higher-order multipole 
form factors. The normalization of all the form factors are unknown except for the charge form factor of G,(0) = 1. 
It should be mentioned that the discussion of all the seven form factors is beyond the scope of this paper, and it will 
appear elsewhere. 


In analogy to the spin-1/2 nucleon and spin-1 deuteron cases, we assume that the charge distribution of the spin-3 
particle, d* (2380), is also directly relate to the matrix element of J? in the Breit frame. For the six-quark system, we 
just consider the quark-quark-photon current 


6 6 
P =X ear a= > 3. (11) 
i=1 i=l 


Thus, we may calculate the matrix element of J” and determine the charge distributions of the d* (2380) in the form 
of 


3 
NO <pyma: |J |p me >. (12) 


Mar =—3 


1 


d* 25 — 
CEQ) = = 


II. CALCULATIONS OF THE d* CHARGE DISTRIBUTION IN TWO SCENARIOS 
A. Scenario A: Hexaquark dominant structure 


Here, we only concentrate on the charge distribution of d*(2380). As mentioned in Refs. ie, our model wave 
function for d*(2380) is obtained by dynamically solving the bound-state RGM equation of the six quark system in 
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the framework of the extended chiral SU(3) quark model, and then successively projecting the solution onto the inner 
cluster wave functions of the AA and CC channels. The resultant wave function of d* can finally be abbreviated to 
a form of 


|W ax(2380)) = a| AA )(sn=(30) + Bl CC )(s1n= 30) (13) 
| [ da (61, &) palé 65) xaa (B) Gan + gelê E) gc lër £) xec(R)) ce | s=), 


where a and @ are the fractions of the AA and CC components in d*(2380), da and ¢c denote the inner cluster 
wave functions of A and C (color-octet particle) in the coordinate space, yaa and ycc represent the channel wave 
functions in the AA and CC channels (in the single AA channel case, the CC component is absent), and ¢aa and 
Cec stand for the spin-isospin wave functions in the hadronic degrees of freedom in the corresponding channels, 
respectively Es. It should be specially mentioned that in such a d* wave function, two channel wave functions are 
orthogonal to each other and contain all the totally anti-symmetrization effects implicitly : 


Unlike in calculations of the decay processes of d* —> daz, d* > NNan, and d* + NNza where the CC component 
does not contribute to the widths, here in the calculation of the charge distribution of the d* both the CC and AA 
components contribute. Considering that A and C are antisymmetric then the charge distribution is 


6 
GQ) = <a (P')| Soi |a(P) > (14) 
i=] 
= 3<d*(P’) | (e345 + e636) | d* (P) >, 
where the superscribe ” A” stands for the scenario A, P' — P = g, Q? = F?, and €3.6 = ł + Be Then 
GBQ’) =3 PUS + IB )OrOy, + PULSE +16 Oc yc | (15) 


where Oa c denote the overlaps of the wave functions of the 3-rd quark (or 6-th quark) which bombarded by photon 
in A and C, and Oy, and Oy, represent the contributions from the AA and CC channel wave functions, respectively. 
ey 3 can be calculated by 


Is = arh AA | ess | AA Yar) (16) 
ISe (1412) CC | e3,6 | CC )iz,z.)- 


Finally, one obtains 


G4’) = e TENN E “Joy, 


- (17) 


where bac are the size-parameters of the A and C systems. As has been discussed explicitly in Ref. ie, the channel 
wave function in the AA channel is written in terms of Gaussian-like wave functions as 


2 


x(r)aa = 3 Feo ( im on) (18) 


where cm and bm can be determined by fitting the channel wave function in this form to our projected model wave 
function calculated before. The normalization condition of the wave function is f d?r | x(r) |?= 1. Thus 


= š b3 b3 27 3/2 b202 
2\— > E men = men -2 
Oxa (Q*) =< xaa lk + G/2) | xaa (k) >= a Emen- (ae ra z) exp| 302 + oe)! |; (19) 


For the CC component, the channel wave function is dominated by the single S-wave Gaussian function 


xoolr) = (tv7) exp |- r2/20]. (20) 


Similarly, the contribution from the CC component can be calculated by 


Oxe(@) =< xoc lË + 4/2) | xoc Ë >= exp | - Pe. (21) 
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B. Scenario B: Dior structure 


We can also calculate the wave function of Diz with (I=1, S=2) by using our chiral SU(3) constituent quark model. 
The obtained mass is about My + Ma — e€ ~ 2175— e MeV (e is binding energy). It may be very close to the threshold 
of NA. Similar to eq. (18), the relative wave function between N and A can also be expressed as 


ius = 


m=1 


2 


Fee (- 3): (2) 


where cl, and b/,, can be determined by fitting the relative wave function in this form to the resultant wave function 
of the Dı in our model calculation. 


We first calculate the charge distribution of Dı2 by assuming it as a NA 6-quark system. The procedure is the 
same as that for the AA system shown in Sec. IIIA. The obtained charge distribution can be written in the following 
form: 


GEQ’; m) = 3] Svs" + Saé], (23) 


where m; stands for the third component of the isospin of Di2 and 


a 2 b3; a 
Sy, a(Q®) = exp | - >] x Ow, (24) 
with 
z b3b3/ 2m 53/2 b22 
Oyna (Q?) =< xna lk + G/2) | xna (k ao Cn a a) ezp - g 4 oa! | (25) 
and 
a 1⁄1 . 3 
E3 t= ( 1, m| €3 | 1, m ) = (z 1, z (mi = 1,0, —1) 
EM 1/7 3 
é™ = (1,m| e6|1,m) = =(5 i -5), (mı = 1,0, —1). (26) 


Now, we calculate the charge distribution of d* which is assumed to have a Dər structure. What we would like to 
see is that if such type of structure of d* has a distinguishable charge distribution compared with the one shown in 
Sec. IIIA. Since the (J)? values of D12 are 1(2)*, the relative motion between D12 and 7 has to be at least P-wave, 
and Djo7 in the isospin space should be decomposed as 


1 
| d >= — [Dit = 177 — Dio(I, = 0)n° + Diall; = —1)rt |. (27) 
V3 
The Jacobi momenta of D;2—7 system are 
a 5 > MrDp,. — Mp Dr 3 = 
P= Bet on iT Mpa T Poa — (28) 


where pr and p,,, are the momenta of t and D12, and q stands for the relative momentum between the two systems. 


From the above equation, one sees that the bombarding effect of photon on the relative momentum q is much 
smaller in the case where Dj is stricken than that in the case where m is hit. This is because of the factors of 


a= mr/(Moyy + mr) ~ 6/100 and b=1—a= Mo, /(Mors + mr) ~ 94/100. 


Although the relative wave function between D12 and 7 is not strictly solved, one can still qualitatively see a general 
character of such a structure by assuming it being a P—wave function 


7 V2b5/? oe 
Xp (P m) = Sqr Vim (Oa) exp(—b*4 2)/2 (29) 
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where Vı,m, (Qg) is a so-called solid harmonics with (l, m) = (1, mz), and b denotes the size-parameter for the relative 
motion between the Dı and 7. The size-parameter b can be adopted in a large range, for instance from 0.6 ~ 6fm, 
and the real relative p—wave function could be a combination of such a function with different size-parameters. Then 
the contribution of the relative wave function between Dj 2 and 7 can be calculated by 


OB (Gm) = (x5 x(k +, m) | Xo, (k, m) ), (30) 


where the subscript Di2 and 7 in Or „ denote the situations when a photon hits Di2 and 7, respectively. The 
phenomenological monopole phranbehrnatiot of the pion charge distribution can be borrowed from Ref. (23) 


5 1 
Sx = FAQ") = T+ (Q2/A2)’ AP = 0.462 GeV?. (31) 


Then, the charge distribution of d* in the Dj27 scenario can be written as 


CB mg lQ) = EJE CLM, a im, (08° (OR (aq, mu) + rO — a)a,m))]. (32) 


Mt, mı 


It should be stressed that the pion contributions from the first and the third terms of eq. (27) canceled each other 
and the one from the second term vanishes since it relates to 7°. Therefore, the charge distribution of d* in this 
scenario only comes from the contributions by D12 and the relative motion between Dım. Averaging over the initial 


states with various magnetic quantum numbers of d*, one finally obtain the charged distribution of d* 


1 
GEQ?) = 5 fois?) x Oxy [exp (— DRT ?/6) + exp (— 0472/6) , (33) 
where 
| 272 272 P 
foal’) = (1 = <7) exp | — =! j: (34) 


C. Numerical results in the two scenarios 


In our calculations, we take by = ba = 0.5 fm, be = 0.45 fm, b; = 0.45 fm, and b=1.2 fm as inputs. The 
probabilities in eq. (13) are a? ~ 0.31 and 8? ~ 0.69. In Fig. 1 we plot the channel wave functions of the AA in both 
the single AA channel (Scenario A1) and coupled AA + CC channel (Scenario A2) approximations, respectively. In 
terms of the same chiral SU(3) constituent quark model, the wave functions of the NA system with (J, S) = (1,2) can 
also been obtained by performing a bound state RGM calculation with a set of slightly varied model parameters whose 
values are still in a reasonable region. It is shown that in such a wave function, the 5S2 component dominates the D12 
state with a fraction of about 94.47%, when the NA system is weakly bound with a binding energy of e = 4.17 MeV. 
The wave function of Dı2 is displayed in Fig. 1. In terms of the wave functions in the Al, A2, and Dj, cases, the 
root-mean-square radii (rms) of the d* can be straightforward calculated. The obtained rms are listed in Table [} 


TABLE I: The calculated rms for d* (2380) in the scenarios Al and A2 and for our D12 (in units of fm). 


Cases d* (2380) Die 
Al A2 


ae | 
Gm) 739 | 


In Fig. 1 and Table |I| one should notice that the AA wave function in the coupled channel approximation is 
normalized to 0.31. Clearly, we see that the wave function of Dı2 is much more extended in radius in comparison 
with that of AA does. This is reasonable, because the energy level of this state is fairly close to the NA threshold, 
the system is almost broken up, namely N and A are ”almost free” and the separation between them becomes rather 
large. The smaller the binding energy is, the larger the size of the system would be. The root-mean-square radius of 
this component is 2.39 fm. This value is consistent with the prediction from Heisenberg uncertainty bil. 
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Scenario A 


roem A 


4 
r(fm) 


FIG. 1: Wave functions. The pink dashed-dotted curve and the black solid curve represent the AA components in the d* wave 
function in the framework of our chiral SU(3) constituent quark model with the single AA (scenario A1) channel only and the 
coupled AA + CC channel (scenario A2), respectively. The blue double-dotted-dashed curve shows the relative wave function 
of the °S2 component in D12 with binding energy of € = 4.17 MeV. 


The charge distributions of the d* can be calculated by Eq. in scenario A. The parameters Cm and bm in 
Eq. (8) can be found in Refs. ud uj. The obtained charge distributions from AA and CC components are 
demonstrated in Fig. 2. The charge distribution of the single AA model (scenario A1) is plotted in Fig. 2 by a pink 
dotted-dashed curve. The black solid and the red dashed curves denote the contributions from the AA and CC 
components (scenario A2), respectively, and the black dotted curve represents the total contribution by summing 
over former two curves. These curves tell us that the contribution from CC component is larger than that from 
the AA component, especially, in the larger momentum transfer region, the contribution is dominated by the CC 
component. It implies that the quark contents in the CC component tend to concentrate in a more compact region 
than that in the AA component. This physics picture coincides with the radii of two components calculated in 
our previous paper (15). Moreover, the curvature of the distribution curve of Al is larger than that in the coupled 
channel case. It indicates that quarks here distribute in a larger region than that in the coupled channel case. The 
size information of d* can also be seen from the slope of the distribution curve at the origin, because such a slope is 
closely related to the radius of the system. The larger the slope at origin is, the larger the radius of the system would 
be. Comparing these slopes with the wave functions shown in Fig. 1, we find that the obtained slopes at the origin 
here coincide with the radii shown in Fig. 1. One sees that the radius resulted from the single channel wave function 
(A1) is larger than one from the two coupled-channel approximation (A2). Moreover, the radius contributed by the 
CC component is smaller than the one from AA component. 


In order to see the size character of d* with a Di2 — m structure, the charge distribution of d* in the scenario B 
is calculated by Eq. (33). The obtained charge distribution curve, depicted by the blue double-dotted-dashed curve, 
is also drawn in Fig. 2 as well. It should be specially mentioned that in our numerical calculation, we do not solve 
the bound state problem for the D,27 system, explicitly. However, since the requirement of the conservations of the 
total spin and parity, the relative motion between Dı2 and m must be at least a P-wave. Therefore, if we ignore 
the component with higher partial wave, which will be greatly suppressed, the true relative wave function would 
be a superposition of the P-wave wave functions with different size-parameters b. Thus, we calculate the charge 
distribution curves with a b value from 0.6 fm to 6 fm. The result shows that those curves almost overlap with 
each other, namely this curve is almost insensitive to the size-parameter b. This is because that the incoming photon 
is absorbed by the Dı2 system, and the induced change of the relative momentum between the Di2 and m is very 


small due to the factor of m,/ (u DaF mr) ~ 6/100. Comparing the curve of scenario B with others of scenario A, 
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Scenario A 
Single channel 
<=. > i = AK 


Coupled-channel 


Q'(GeV’) 


FIG. 2: The obtained charge distributions of d*. The solid(black), red(dashed), and black(dotted) curves stand for | AAxa >, 
| ccxe >, and the total contributions in the coupled channel approximation, respectively. The pink dotted-dashed curve 
represents the single channel case of scenario A. The blue double-dotted-dashed curves stand for the results of Dj27 scenario 
with e ~ 4.17 MeV. 


we see that the curves for the Dior scenario decrease and go to zero much faster than those for the scenario A. A 
much larger slope at the origin means that the radius of the Dior system is much larger in comparison with those in 
scenario A. 


From our numerical calculation, we find that the ratios of the slopes of the curves at the origin in scenarios A1, 
A2, and B are 


Al 2 A2 2 2 
R= e ) /(- 2E ) /(- 29 ) / 


It should be mentioned that the contributions to 2e, (2) from AA and CC components in the coupled channel 
approximation of Scenario A2 are 0.64 and 0.81, respectively. The above obtained ratios tell us that the slop of the 
charge distribution of the d* in scenario B at Q? = 0 is about 4.8 times larger than the one in Scenario A2, and 
the corresponding charge radius of scenario B is about 2.2 times larger than the one in scenario A2. This feature is 
compatable with that discussed in Ref. RA. Finally, the very sharp charge distributions (or large charge radii) of 
d* (2380) in the scenario B is mainly dominated by the very broad wave function of the obtained D12. Therefore, we 
conclude that the two scenarios, A2 and B for the d* (2380) give very different descriptions for its charge distribution 
and its charge radius. 


= 2.30 : 1.45 : 6.86. (35) 
Q?=0 


IV. SUMMARY 


We have calculated the charge distribution of the d* in the two scenarios; one is a hexaquark dominant picture and 
another is a Dj27 resonant picture. In the first picture, we show the total charge distribution and both contributions 
from its AA and CC components. In order to make a comparison, the result of the single AA channel is also 
shown. Comparing the predictions of the two scenarios, we see that the charge distribution from the Dj27 system is 
remarkably different from the scenario A2, and consequentially, the charge radius of the the scenario B is obviously 
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larger than that of the scenario A2. 


We now expect a series of experiments which may be able to test different interpretations of d* in future. Although 
the direct ed* scattering measurement may be hard to carry out, one may consider the d* form factors in the time-like 
region. For example, the production of the final d*d* pair in the e+e~ and pp annihilation processes. It is our hope 
that the future upgraded BEPC, Belle and Babar and experiment at Panda with high luminosity may provide a test 
for different theoretical understandings. 


Acknowledgments 


We would like to thank Heinz Clement, Qiang Zhao, and Qi-Fang Lü for their useful and constructive discussions. 
This work is supported by the National Natural Sciences Foundations of China under the grant Nos. 11475192, 
11475181, 11521505, 11565007, and 11635009, and by the fund provided to the Sino-German CRC 110 “Symmetries 
and the Emergence of Structure in QCD” project by NSFC under the grant No.11621131001, the IHEP Innovation 
Fund under the grant No. Y4545190Y2. F. Huang is grateful for the support of the Youth Innovation Promotion 
Association of CAS under the grant No. 2015358. 


any 


M. Bashkanov et al., Phys. Rev. Lett. 102 052301 (2009). 

2] P. Adlarson et al., Phys. Rev. Lett. 106, 242302 (2011); P. Adlarson et al., Phys. Lett. B 721, 229 (2013); P. Adlarson et 
al., Phys. Rev. Lett. 112, 202301 (2014). 

3] A. Abashian, N. E. Booth, and K. W. Crowe, Phys. Rev. Lett. 5, 258 (1960); N. E. Booth, A. Abashian, and K. M. Crowe, 
Phys. Rev. Lett. 7, 35 (1961); F. Plouin et al., Nucl. Phys. A 302, 413 (1978), J. Banaigs et al., Nucl. Phys. B 105, 52 
(1976). 

4] M. Bashkanov, H. Clement, T. Skorodko, Eur. Phys. J. A51, 87 (2015). 

5] COSY confirms existence of siz-quark states, CERN COURIER, Vol. 54, No. 6, p6, July 23, 2014. 

6] H. X. Chen, W. Chen, X. Liu and S. L. Zhu, Phys. Rept. 639, 1 (2016). 

7] J. Dyson, Phys Rev. Lett. 13, 815 (1964). 

8] A. W. Thomas, J. Phys. G 9, 1159 (1983). 

9] M. Oka and K.Yazaki, Phys. Lett. B 90, 41 (1980). 

10] T. Goldman et al., Phys. Rev. C 39, 1889 (1989). 

11] X. Q. Yuan et al., Phys. Rev. C 60, 045203 (1999). 

12] M. N. Platonova and V. I. Kukulin, Phys. Rev. C 87, 025202 (2013); Nucl. Phys. A 946, 117 (2016). 

13] H. X. Huang et al. Phys. ReV. C 79, 024001 (2009); Phys. ReV. C 89, 034001 (2014). 

14] A. Gal and H. Garcilazo, Phys Rev. Lett. 111, 172301 (2013); Nucl. Phys. A 928, 73 (2014); A. Gal, Acta Physica Polonica 
B 47, 471 (2016). 

15] Fei Huang, Zongye Zhang, Pengnian Shen, and Wenling Wang, Chin. Phys. C39, 071001 (2015), Fei Huang, Pengnian 
Shen, Yubing Dong, and Zongye Zhang, Sci. China 59, 622002 (2016), and references therein. 

16] Yubing Dong, Pengnian Shen, Fei Huang, and Zongye Zhang, Phys. Rev. C91, 064002 (2015). 

17] Yubing Dong, Fei Huang, Pengnian Shen, and Zongye Zhang, Phys. Rev. C94, 014003 (2016). 

18] Yubing Dong, Fei Huang, Pengnian Shen, and Zongye Zhang, ”Decay width of d*(2380) = NNvz process in a chiral 
constituent quark model”, arXiv: 1702.03658v2 [nucl-th]. 

19] The WASA-at-COSY Collaboration, ”Isoscalar Single-Pion Production in the Region of Roper and d*(2380) Resonances” 
arXiv: 1702.07212V1 [nucl-ex]. 

20] M. Bashkanov, Stanley J. Brodsky, and H. Clement, Phys. Lett. B 727, 438 (2013). 

21| H. Clement, ”On the History of Dibaryons and their Final Observation”, Progress in Particle and Nuclear Physics, 93, 
195 (2017); arXiv:1610.0559v1 [Nucl-ex]. 

22] N. Barik, S. N. Jena, and D. P. Rath, Phys. Rev. D41, 1568 (1990). 

23] S. R. Amendolia et al., Nucl. Phys. B277 (1986) 168. 

24] A. Gal, ” The d* (2380) dibaryon resonance width and decay branching ratios”, arXiv:1612.05092, Phys. Lett. B, In Press. 


